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Fran^oise  Cornu  ,  Bernard  Jancovici  ,  and  Lesser  Blum 

Abstract 

Using  a  new  method,  we  find  that  the  two-dimensional  one-component 
plasma  in  an  inhomogeneous  background  is  a  solvable  model,  in  equilibrium 
classical  statistical  mechanics,  for  the  special  value  f  »  2  of  the  coupling 
constant,  for  a  larger  class  of  background  shapes  :  the  n-body  densities  con 
be  explicitly  computed.  In  particular,  we  can  deal  with  a  doubly  periodic 
background;  this  is  a  classical  model  for  a  crystal  made  of  fixed  ions  and 
mobile  electrons.  At  F  -  2,  this  system  is  conductor ;  the  correlations  have 
a  fast  decay,  and  the  Stillinger-Lovett  screening  sum  rule  Is  obeyed. 

KEY  WORDS  ;  One-component  Inhomogeneous  plasma;  doubly  periodic 
background;  conducting  phase. 
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t.  INTRODUCTION 

In  statistical  mechanics,  it  is  obviously  of  interest  to  have  exactly 
solvable  models  for  Coulomb  systems  (plasmas,  electrolytes,  metals  etc  ) 
The  simplest  model  of  a  Coulomb  system  is  the  one-component  plasma 
(jellium)  :  identical  charged  particles  move  in  a  rigid  charged  background 
which  ensures  overall  neutrality.  In  two  dimensions,  the  Coulomb  potential 
between  two  particles  of  charge  e  at  a  distance  r  from  one  another  Is  -e^ 
in(r/\),  where  \  is  an  arbitrary  length  scale,  and  the  dimensionless 
coupling  constant  is  r  «  §e^,  where  ^  is  the  inverse  temperature;  for  the 
special  value  F  «  2,  it  has  been  previously  found  that  the  equilibrium 
classical  statistical  mechanics  of  the  two-dimensional  one-component 
plasma  can  be  worked  out  exactly  for  several  kinds  of  background  charge 
distributions  :  one  is  able  to  obtain  the  n-particle  densities.  Besides  the 
simplest  case  of  a  uniform  backgroura/*'^^  essentially  one  could  deal  with 
a  background  charge  density  depending  on  one  space  coordinate^^^;  this 
covers  a  variety  of  charged  interfaces  (electrical  double  layers)  of  interest 
to  electrochemists. 

This  previous  work  left  unsolved  the  important  case  of  a  doubly 
periodic  background,  in  the  present  paper,  we  solve  this  case,  using  a  more 

(4) 

general  new  method;  a  preliminary  account  has  been  given  by  two  of  us  . 
Thus  we  have  an  explicit  solution  for  a  model  which  can  be  understood  os 
made  of  mobile  (classical)  ‘electrons’  interacting  between  themselves  and 
with  a  lattice  of  extended  fixed  ‘Ions*;  this  one-component  plasma  in  a 
periodic  background  con  also  be  regarded  as  a  two-component  plasma  in 
which  the  particles  of  one  species  have  been  fixed  on  a  lattice.  Like  the 
symmetric  two-component  plasma,  the  present  model  is  expected  to 
undergo  a  Kosterlitz-Thouless  phase  transition  between  a  low-temperature 
dielectric  phase  and  a  high-temperature  conducting  phase,  and  this 
transition  is  actually  seen  in  computer  simulotions^^'^^  Here,  we  show 


uaiBuauiBvu 
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that,  in  our  system,  cntena  which  characterize  a  conductor  are  satisfied 
the  correlations  have  a  fast  decay  at  large  separations,  and  the 
Stillinger-Lovett  sum  rule^^^  is  obeyed.  According  to  these  criteria,  at  f  ■ 
2,  the  system  is  in  its  conducting  phase 

The  paper  is  organized  as  follows.  In  Section  2,  the  general  method  is 
reviewed.  In  Section  3,  it  is  shown  that  this  method  provides  a  simpler 
approach  to  the  known  case  where  the  background  density  is  inhomogeneous 
in  one  direction  only.  The  doubly  periodic  background  is  discussed  in  Section 
4 :  we  compute  the  n-particle  densities  and  discuss  sum  rules. 

2.  HETHOO 

2.1.  n-particle  densities 

We  start  with  N  particles  of  charge  -e  in  some  background.  The 
position  of  the  ith  particle  is  -  (x^,  y^);  we  shall  also  use  the  complex 

number  2^  •  x^  •»  iy^.  The  Hamiltonian  is 
N 

H-e^£  V(rt)-e^  £  indz^  -  zjIA)  (2.1) 

i-1  l<i<)<N 

where  e^V  is  the  background-particle  interaction,  and  therefore,  for  an 
inverse  temperature  g  such  that  r  >  ge^  -  2,  the  Boltzmann  factor  is 

exp(-9H)  -  C  Idet  {exp  (-V(rj)l  zj  )i  j- 1 , .  ,  ^2  2) 


where  C  is  a  constant. 

In  the  simple  case  of  a  background  potential  of  circular  symmetry,  V(r) 
•  V(r),  the  functions  expC-Vfr)]  z^~  *  ere  mutually  orthogonal,  and  (2.2)  has 
the  same  form  as  the  squared  wavef unction  of  a  system  of  independ|nt 
fermions;  to  dee)  with  the  Slater  determinant  which  appears  in  (2.2)  is  a 

*  ^  na-A-— V  y,  ^7.  -  2  ^  ^  j ' — A_ 
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standard  problem,  and  it  is  easy  to  compute  the  n-particle  densities  All  the 

zxjfiK — Crutsjii't"  u«  i  7 

previously  solved  cases  could  be  obtained  by  starting  with  a  circular 
geometry  and  taking  a  suitable  limit. 

In  the  present  paper,  we  want  to  consider  more  general  forms  of  the 
background  potential,  and  the  functions  exp(-v(r)]  ore  not  necessarily 
mutually  orthogonal.  However,  we  con  choose  on  orthogonal  basis  fj(r)  for 

the  space  of  these  functions,  and  rewrite  (2.2)  as 


exp(-§H)  -  C  ldet{V^(r^)}^  ^  ^ 


since  the  new  determinant  is  proportional  to  the  former  one.  It  is  then  easy 
to  show  that  the  n-particle  truncated  densities  can  be  expressed  in  terms 
of  the  projector 


<r,IPIr2>-£ 


Vj(r,)Vj(r2) 

I  d  r1¥j(r)l^ 


•I 

1 


p(r)  -  <rlP|r> 

9®V,,  Tj)  -  -l<  r,  IPI  Tj  >1* 

r, . r„)-  -I  <r,  IPIr,  >  ..<r|  IPIr,  > 

(I,  V-V  ’  »  "  ' 


Where  the  summation  runs  overall  cycles  (I,  t2-ln)  hullt  with  {1.2.  ...n}.  in 

the  thermodynamic  limit,  the  functions  exp(-V(r))z^”  *  span  the  subspace  of  rU-e. 
Hilbert  space  defined  by  the  entire  functions  of  z  -  x « iy  times  exp(-V(r)I 
and  P  becomes  tN  projector  on  that  subspace  (of  course,  this  is  an  intrinsic 
definition  of  P,  indeperxtont  of  the  choice  of  the  orthogonal  basis  Thus, 


the  problem  of  obtaining  the  n-particle  densities  is  reduced  to  computing 
the  projector  P 

In  the  simplest  case  of  a  uniform  background  density  the 
background  potential  V(r)  can  be  chosen  as  1/2  it^r^  (plus  some  irrelevant 
constant),  H'j  -  exp(-V)z^~*  and 


<  TjlPlfj  >  *  Po 

Q(r)  -  Po 

•’2^  •  -Po  exp(-npolr,  -  r/) 


In  the  genera)  case,  the  background  density  Pg(r)  can  be  considered  as 

being  the  sum  of  a  uniform  contribution  plus  a  non  uniform  modulation 
p(r).  Correspondingly,  the  background  potential  V(r)  can  be  chosen  of  the 
form  VjjCr)  ♦  0(r),  where  V^Cr)  •  (1/2)  and  Ae^Or)  •  2rtp(r).  As  a  first 
step  towards  the  computation  of  the  projector  P  on  the  space  of  the 
functions  exp(-fif(r)  -  Vp(r))z^,  J  €  N,  it  will  turn  out  to  be  convenient  to 
replace  the  by  another  basis  for  the  entire  functions  Pj^fz)  - 

exp{-(1/2)np^z-(k/Hpg)f},  k  €  R.  The  are  indeed  such  a  basis,  since 
00 

z”  -  n  ’^2“"  I  dt  Hn(t)  exp(-(z  -  X^]  (2.7) 


where  the  H^(t)  ere  Hermite  polynomials,  and  (2.7)  becomes  a  superposition 
of  functions  through  a  rescaling  of  z  and  t.  The  basis  exrt-Vo(r))z^  is  then 
replaced  by 

8xrf-VQ(r)Jp.  (z)  -  exp(-iHp^xy)  exp(-k^/4ng^)  exp{-n9^x-(k/2j(gQ)f  ♦  iky} 
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Actually,  since  (2,8)  will  be  used  for  defining  the  projector  P  and  thereafter 
computing  the  densities  (2.5),  we  can  omit  in  (2.8)  the  normalization  factor 
exp(-k^/4npg)  (this  leaves  the  projector  unchanged)  and  the  phase  factor 
exp(-ii(9Qxy)  (this  leaves  the  densities  unchanged  because  of  their  cyclic 
structure). 

Therefore,  an  alternative  definition  of  P  is  to  take  it  as  the  projector 
on  the  space  of  the  functions 

f  k(r)  -  exp(-0(r)]  expf-rtgufx  -  -!^  f  ♦  iky]  (2.9) 

2«go 


For  the  potentials  0(r)  which  will  be  considered  here,  strictly 
speaking  the  functions  (2.9)  do  not  belong  to  Hilbert  s  space,  because 
lexp(iky)i  does  not  decrease  at  infinity.  However,  these  functions  do  form  a 
basis  in  the  sense  of  distributions,  just  like  the  plane  waves  in  quantum 
mechanics. 


Let  us  remark  that  the  background  potential  V(r}  is  not  uniquely 
determined  by  the  background  density  p^Cr).  In  the  thermodynamic  limit. 


v(r)  keeps  a  memory  of  the  boundary  conditions  even  after  these  boundaries 
have  receded  to  infinity.  Die  n-particle  densities  however  should  depend 
only  on  Q^(r)  for  a  system  with  screening  properties  which  prevent  the  bulk 


from  being  affected  by  infinitely  remote  charged  boundaries,  it  is 
satisfactory  to  check  explicitly  this  independence  upon  the  choice  of  V(r), 


to  which  it  should  be  always  possible  to  add  an  arhitrary  harmonic  function. 
Actually,  since  the  confinement  of  the  particles  must  be  preserved,  the 
total  background  potential  must  increase  fast  enough  at  infinity,  and  we 
shall  only  consider  the  addition  to  V(r)  of  a  harmonic  function  of  the  form 


f(r)  =  a2(x^-y^)  ♦  h^xy  +  a,x  ♦  b,y  ♦  c,  with  coefficients  and 
sufficiently  small  absolute  value,  furthermore,  the  term  b^xy  can  be 
removed  by  a  rotation  of  the  axes.  The  potential  f(r)  can  be  interpreted  as 
determined  by  suitable  external  electrodes. 

It  is  then  easy  to  see  that  ^(z)  can  be  chosen  again  of  the  form 

exp{-a{z  -  (2ar'(k  *  x  *  where  a,  y,  6  are  real  constants,  and  by  a 
suitable  choice  of  these  constants,  f(r)  can  be  cancelled,  except  for 
irrelevant  normalization  and  phase  factors.  Therefore,  the  projector  P  can 
be  left  unchanged. 

Incidentally,  the  basts  (2.9)  cah  be  directly  obtained  by  choosing  the 
potential  contribution  from  the  uniform  background  as  and  taking 

exp(k2)  as  the  basis  for  the  entire  functions.  Such  a  will  be  obtained  if 
we  reach  the  infinite  system  limit  starting  from  a  strip  geometry  rather 
than  from  a  circular  one. 

2.3.  naonetlc  analoou 

The  arbitrariness  in  the  choice  of  our  background  potential  V(r)  has  a 
quantum-mechanical  analog  ;  a  gauge  transformation  in  a  magnetic  problem. 

Let  us  consider  a  particle  of  mass  m  and  charge  q  moving  in  the  xy  plane  and 

(B) 

submitted  to  a  uniform  magnetic  field  B  parallel  to  the  z  axis  .  The  ground 
state  is  infinitely  degenerate.  In  the  gauge  where  the  vector  potential  is  A 
<-  (1/2)6  X  r  ,  a  basis  for  the  ground-state  wavefunctions  is 

ex|](-(qe/4}r^);^,  j  c  N;  in  the  gauge  where  the  vector  potential  is  A  «  j  Bx 
(J  is  the  unit  vector  along  the  y-axis),  a  basis  for  the  ground-state  wave 
functions  is  exp{-(qB/2)  [x  -  (k/qD)f  *  iky),  k  €  R.  Obviously,  our  change 
from  the  basis  exp(-(  I /2)it^r^)z^  to  the  basis  exp{-rt^x  -  (k/2icgg)f  *  iky) 
is  exactly  of  the  same  form.  If  a  Slater  determinant  Is  built  with  the 
ground-state  wave  functions,  the  corresponding  n-body  densities  must  be 
gauge-independent,  just  as  our  n-body  densities  ore  independent  of  our 
choice  of  basis. 


0 


3.  BACKGROUND  INHOnOGENEOUS  IN  ONE  DIRECTION. 


We  now  revisit  the  case  of  a  background  density  depending  on  one 


coordinate  only^  pg(x).  The  potential  9i  can  be  chosen  as  0(x),  and  we 


have  at  hand  the  orthogonal  basis  ;  (2.9)  is  orthogonal  because  of  the 
plane-wave  factor  exp(lky).  Adapting  (2.4)  to  the  case  of  a  continuous  index 
k,  with  the  usual  prescription  that  the  summation  becomes  (L/2n)  |  dk  as 
the  length  L  of  the  system  in  the  y  direction  becomes  Infinite,  we  find 


<  r,  IPI  Tj  >  -  ex()(-0(x,)  - 


ilk 

•00  2n 


exp(1k(y,-y2)l  expl-rtg^(x, - k— )^  ♦  (Xj  -  — k— )^]} 

2n^ 


f®®  dx  exp(-20(x)  -  2ne^(x  -  -k_)^] 

-o®  2n^ 


Using  (3.1)  In  (2.5),  we  retrieve  st  once  the  ueiiiltlic  of  Ref .3. 


Some  special  care  must  be  exercised  for  dealing  with  the  case  where 
the  particles  are  confined  to  the  half  space  x  >  0  by  an  Impenetrable  wall  at 
X  «  0.  Then.  In  (3.1).  the  range  of  x  must  be  restricted  to  x  >  0.  and  the 
range  of  k  must  be  restricted  to  K  >  0.  This  Is  shown  as  follows.  Since  we 


have  already  taken  the  limit  of  on  Infinite  system,  we  shall  start^  ^  with  a 
system  In  which  on  Impenetrable  barrier  occupying  the  region  -  I  <  x  <  0 
separates  the  plasma  Into  two  regions  x  <  -  I  and  x  >  0.  Ttie  impenetrable 
wall  system  will  be  obtained  by  taking  the  limit  I  oo  in  such  a  way  that 


the  remote  regions  x>Oandx<-lno  longer  see  each  other.  The  reason  for 

u  O 

Which  the  values  k  <  0  ore  suppressed  In  (3.1)  will  bo  aeon  to  be  that  the 


norm  In  the  denominator  has  a  contribution  from  the  remote  region  x  <  -  I 
which  becomes  Infinite  for  k  <  0,  In  the  limit  I  -►  oo;fti  the  contrary,  this 


contribution  vanishes  for  k  >  0,  In  the  limit  i  -►  ».  The  Independence  of 


the  regions  x  >  0  end  x  <  -I  oan-M  achieved  by  requiring  that  each  of  them 


potential  symmelncj^  with  respect  to  the  banner,  V(x)  =  ♦  0(x)  for 

X  >  0,  V(x)  =  +  oo  for  -i  <  X  <  0,  V(x)  =  V(-x  -i)  for  x  <  -i.  Then, 
remembemng  that  0(x)  is  defined  in  every  region  as  V(x)  -  x^,  and 

changing  -x  -i  into  x,  we  can  rearrange  the  denominator  of  (3.1)  as 


foo 

dx  exq(-20(x)  -  2n9n(x  ♦  — )^] 

0  2ngo 


♦  I  dx  exp(-2£f(x)  -  2itgn(x - ^ — )^] 

0  2ngo 


(32) 


In  the  limit  i  the  first  term  of  f(k)  (which  is  the  contribution  from 
the  region  x  <  -i)  diverges  if  k  <  0,  vanishes  if  k  >  0  As  a  consequence, 
(3. 1 )  becomes 


<  r,  IP|  Tj  >  -  exp(-0(x,)  -  0{x^)] 


X 


f 

OO 

^0  2rt 


exp(ik(y,-y2)]  exp{-«9j(x,  -  ♦  (x^  -  — !^)^]} 

_ _ 2"^ 

f“  dx  exp(-20(x)  -  2K9o(x  - 

0  2j19o 


,x,.X2>0 

(33) 


(2) 

(3.3)  is  In  agreement  with  the  results  about  a  hard  wall  carrying  a 
surface  charge  density  etr,  in  which  case  0M  »  2JKrx.  x  >  0;  0(x)  -  0,  x  <  0 


4.  DOUBLY  PERIODIC  BACK6R0UND 


In  this  Section,  which  is  the  core  of  the  present  paper,  we  study  the 
cose  where  the  background  density  Is  doubly  periodic.  Thus,  we  consider  a 
doubly  periodic  background  potential  modulation  0(r) : 


0(r  ♦  n  a  ♦  m  D)  -  0(r),  n,  m  €  Z 


(4.1) 


.  » •  »  «  •  •  ^  •  •.  -  •  •  •  .  .*  ,  *  •  «  ^  •*  •  ■" 

✓  V  v\-  V  ^  **'■  V*" 

^  1^,  i  iftnTi.  Jlni'  S  A  ^  ^  m.  A  m  A  —  M  ■»  ."T  ^  *  A..*  ^  w  ««.  R  —  -<t  ^  . 


The  unit  cell  is  a  parallelogram  built  with  the  vectors  a  and  b,  of  area  ab 
sin  where  ^  is  the  angle  between  a  and  b.  The  system  is  neutral,  with  a 


particle  density  equal  to  the  average  of  the  total  background  density  * 
(2j(r'A0(r).  In  order  to  mimic  a  simple  crystal  of  extended  fixed  ions  'tnd 
mobile  electrons,  of  opposite  charges,  we  take  =  (ab  sin  {)T\  which 


means  there  is  one  particle  per  unit  cell. 


Although  the  functions  (2.9)  now  are  not  orthogonal,  they  are  a  good 
storting  point  for  computing  the  projector  P.  Choosing  the  y  axis  along  the 
period  vector  b,  and  defining  C  €  [0,  I]  and  n  integer  by  k  -  2n(C  *  n)/b,  we 


can  rewrite  as 


Vr  «(r)  -  exp(-flf(r)l  exj<-Kfo(x-  f  *  2ni(C  ♦  n)  -i-  J  (4.2) 

Qo^  b 


As  a  consequence  of  the  perioducity  of  0  along  the  y  axis. 


1“  (H n- <rt 't’c,  n(r) «  9«  -  O 


Furthermore,  if  the  unit  cell  is  a  rectangle  (^  -  n/2),  (Pob)"’  ■  a.  and  the 
periodicity  of  0  along  the  x  axis  ensures  that  depends  on  x  only 

through  x  -  na;  this  suggests  to  introduce  the  Bloch  functions 


ti^r)  -  £  exp(-2niTjn)  .  C,  t)  c  [0,  1] 


which  do  have  the  desired  orthogonality  property 


‘•jiVv  w  v.v..<s 


1  dx  «  S(T1  -  T]')  (4  5) 

The  argument  can  be  easily  extended  to  the  more  general  case  where  the 
unit  cell  is  a  parallelogram,  by  introducing  (dimensionless)  oblique 
coordinates  (X,  Y)  defined  by  r  =  Xa  ♦  Yb.  Multiplying  ^  by  an  irrelevant 

phase  factor  exp{-in(a/b)  cos  4.1^  ♦  (C  ♦  n)^]},  we  obtain 

H'/  n(r)  -  exp(-0(r)]  exrf-  A  (X  -  C  -  n)^  ♦  2«i(C  ♦  n)Y]  (46) 

Where  t  -  (b/a)  exp{1(f  -(rt/2)5;  with  these  functions,  which  have  the 

same  form  as  (4.2),  we  can  proceed  as  above,  in  oblique  coordinates  in  the 
general  case. 

Thus,  the  ¥,  are  orthogonal : 

1  dr 


where 


oo 

f(C,  Ti)  -  J- 1  exp<2j(liiN)  J  dX  ex|<-  X  (X  -  -  JL  (x  -  ;  -  N)^l 

00  N  ^  ^ 


X  1  dY8xp(-20(X,Y))exp(-2niNY) 
0 


These  orthogonal  ¥.  can  be  used  for  building  the  projector  (2.4),  with  the 
s#M 


result 


.1  1 


<  r,  IPl  >  -  exp(-0(r,)  -  ^(fj)]  f  dC  f  dij  — ^ —  Y,  ®xp(2niT|(m  -  n)] 
*  ■‘0  0  f(C,Ti)  nm 
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i 


i 

t 


.1 

A 

A 


( 


3 


H 


X  exp{-  (X,  -  c;  -  n)^  -  i  (x^  -  C  -  m)^  ♦  2ni[(C  >  n)Y,  -  (C  *  m)Y2]}  (4  9) 

A  more  compact  form  can  be  obtained  by  using  the  Poisson  identity 

£  exp(-X  ♦  2j(i2N]  -  </t*  £  ex()(-nt*(z  -  N)^]  (4. 10) 

N  N 

in  (4  8),  and  also  in  (4.9)  where  we  set  m  >  n  «  N.  The  result  is 
T))  -  -L  A*  j“  dX  J“  dY  exp(-2<f(X,  Y)] 

^  -00  -OO 

X  exrt-  S  {%-lf  -  nt*(Y  -  ii)^-2*i(X  -  CKY  -  t\)]  (4.1 1) 

and 

<  r,  IPI  >  -  ex|)(-i?(r,)  -  ^(r^)] 

X  dcl^dti  -J— 8xj)(-  S  (X,-C)*-wt*<Y2-‘n)^-2«1(X2-0(Y2-Ti) 

-OO  -OO  ^ 

♦  2«H(Y, -Yj))  (4.12) 

Thus,  we  have  obtained  an  integral  representation  of  the  projector  P;  when 
used  in  (2.5)  it  gives  the  n-body  densities.  Some  symmetries  of  (4.12)  are 
hidden;  for  instance  <  r,|Pj  >  -  <  Tj  IPI  r,>*.  Other  ones  are  apparent;  for 
instance,  the  function  f(^,  t))  is  doubly  periodical  with  period  1.  and  this 
ensures  that  the  densities  have  the  same  periodicity  properties  as  the 
background. 


As  an  illustration,  we  consider  the  case  of  a  square  unit  cell  (a  -  b,  p 
n/2,  X  -  1)  with  the  simplest  choice 


A 

'lO 

iS 
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exp(-20(r)]  =  t  +  X(cos  2nX  +  cos  2itY)  ,  IX|  ^  1  /2  (4. 1 3) 


The  total  background  density  Is 


9g(r)  =  9q  ♦  (2h)~'a0 


(4.14) 


and  from  (4. 1 1 )  and  (4. 1 2).  we  find  for  the  particle  density 
p(r)  -  <  rIPI  r  > 

-  «pt-2«r)l  foe  J"  1) 

®®  1  ♦  xe'^^'^fcos  2nc  ♦  cos  2mi) 

(4.15) 

The  potential  modulation  0(x,  Y),  the  background  density  pg(x,  Y)  and  the 

particle  density  q(x.  Y)  are  displayed  in  Fig.!,  for  \  •  .49.  With  this  choice 
of  a  large  amplitude  K  the  background  density  is  a  rather  tormented 
lanscape;  the  particle  density  tries  to  follow,  but  it  does  not  quite  succeed 
and  It  exhibits  much  smoother  oscillations. 

Another  representation  of  gtr)  might  be  of  interest.  The  normalization 
factor  (4.6)  can  be  written  as 

f(C,  T|)  -  -L  f'dx  J  dY  t,(X,  Y)P  (4.16) 

Pq  V  0 


and  therefore 


«  0  i;  m,  i;  DY,  Y,)|* 


(4.17) 


If  we  restrict  ourselves  to  the  cose  of  a  rectangular  unit  cell,  t  is  real,  and 


M 


kz  =  exp(-20(X,Y)]  I  exp(-  ?(X-C-n)^-J(X-C-m)^  ♦  2ni(n-m)V-T))] 

'  '  m.n  z  z  T 

(4.18) 

The  sum  on  n  and  m  can  be  replaced  by  a  sum  on  y  =  m  -  n  and  v  =  m  ♦  n, 
with  y  and  v  of  the  same  parity.  The  contributions  from  (y,  v)  even  and  (y,  v) 
odd.  respectively,  can  be  expressed  in  terms  of  Jacobi  theta  functions  The 
result  is 

Kv*’  Wl  (§)'"l8j(2t|  -  2^,  -  1..  |).82(2i|  -  2!^,  2) 

(4.19) 

Where  the  8  functions  are  defined  by 

03(x.  t)  -  £  oxp[-rttn^  ♦  2«inK)  •  exp(-(n/t)(x  -  n)^] 

n  n 

02<>‘*  •  Z  0xp(-«t<n  ♦  2Xi(  n  ♦  -^x) -  £  (- l)"exp(-(n/t)(x  -  n)^] 

04(x,  t)  -  £  exiX-ntn^  ♦  2xin  (x  ♦  ^)]  -  ^  exp(-(n/t)(x  -  n  -  ^)*] 

The  particie  density  is  obtained  by 


tiefn9(4.l9)  in ‘(41 7). 


i  4.2.  Decau  of  the  correlations 

[  The  decay  of  the  truncated  densities  (2.5)  at  large  separations  is 

[  faster  than  any  inverse  power  law.  This  can  be  seen  as  follows.  The  decay  of 

y  the  densities  is  governed  by  the  decay  of  P;  the  0Mpresgl0ir(4.l2)  uf  wmcir 

<  r,  IPI  Tj  >  -  exp(-0(r,)  -  0{r^)]  exp(2j(iX,(Y,-Y2)] 


••  .V 


I 


fOO  f vO 

dC  - - -  expt-Hc;2- nxV  -  2niCTi  ♦  2nij;(Y,-YJ 

-oo  -oo  f(C*X,,Ti*Y2)  ^  ^ 

-  2niTi(X,-X2)] 


(4.21) 


Let  us  study  the  behavior  of  (4.21),  for  a  fixed  value  of  r,,  as  recedes  to 
infinity.  Since  f  is  penodic]|J_in  t)  ♦  Y^,  we  can  define  Y^  by  Y^  =  Y^  +  n,  Y^  e 
[0,  tin  integer,  and  replace  by  Y^  in  f  If  we^e|RfeftC«^elves  to  a 
fixed  value  for  Yj,  i.e.  if  recedes  to  infinity  by  integSr^steps  of  Yj,  the 
integral  in  (4.21)  is  the  Fourier  transform  of  a  well-behaved  function  of  C 
and  t),  and  it  decays  faster  than  any  inverse  power  law. 

The  resulting  fast  decay  of  the  correlations  is  one  of  the  criteria 
which  is  believed  to  characterize  the  conducting  phase. 


f; 


4.3.  Sum  rules 

The  one-  and  two-particle  densities  can  be  shown  to  obey  several  sum 
rules  which  characterize  a  conductor. 

Neutrality 

The  averages,  on  a  unit  cel),  of  the  particle  density  and  of  the  background 
density  are  equal : 


(«( 


dY  9(X,  Y)  -  Po 


(4.22) 


This  sum  rule  is  satisfied  by  (4.17). 

Screening  of  a  particle  of  the  system 

This  screening  rule  means  that  a  particle  of  the  system  induces  a 
polarization  cloud  of  exactly  opposite  charge  : 


V, J"*'  ’."^'  •’^'  ’■'^  U"* 


Uf^  p‘^’(r,,  r^}  -  -9(r,) 


(4.23) 


The  structure  of  (2  5)  ensures  that  this  rule  is  obeyed,  because  of  the 
closure  property 


Jdr2  <  r,|P|  Tj  ><  r2lPI  r,  >  -  <  r,IP|  r,  > 


(424) 


The  rule  can  also  be  checked  explicitly  on  the  representation  (49),  after 
some  algebra. 


This  screening  rule  means  that  an  external  infinitesimal  test  charge 
induces  in  the  system  a  polarization  cloud  of  exactly  opposite  charge. 
Through  linear  response  theory,  this  statement  becomes  a  sum  rule  for  the 
truncated  two-body  density,  the  Camle  and  Chan  generalization^ of  the 

Stllllnger-Lovett  rule^^^  to  an  inhomogeneous  system,  written  here  for  two 
dimensions : 


-pfdr,  fdrj  in  rj  S(r,,  fj)  -  i 


(4.25) 


where  S  is  the  total  charge  structure  factor 


S(r,.  rj)  -  e^  g®’  (r,.  r2)  ♦  g(r,)  6  (r2-  r,)l 


(426) 


In  the  present  cose  of  a  periodic  system.  (4.25)  can  be  written  in  other 
forms.  Because  of  its  periodicity  properties.  S<r,.  r,>  !•  comple^ 
described  by  its  double  Fourier  transform 

B.r,  ♦  I1t.(r2-r,)l  S(r,.  rj  (4  2?) 


where  U  means  that  the  integration  domain  of  r,  t  ■  the  unit  cell,  of  area  A  = 
ab  sin^),  6  IS  a  vector  of  the  reciprocal  lattice  From  (4  23),  one  .deftre^  3 
Sg(6)  =  0  Then,  (425)  can  be  reexpressed  in  terms  of  Sg(k)  and  of  the 

Fourier  transform  2n/k^  of  -in  r  as 


So(k) 

2nB  lim  =  1 
k-0  k^ 


(4  28) 


In  other  words 


I  [  dr,  ldr2  exp(fk.(r2-r,)]S(r,,  r-)  -s.  (4  29) 

Expanding  the  exponential  in  (4.29),  we  find  the  generalization  to  a  periodic 
system  of  the  Stillinger-lovett  second  moment  rule 


1  nr,  foTj  (fj-r.V  (fj-  r,), 


(430) 


is^8(4.: 


where  «,  y  -  1,  2  are  cartesian  compon8nts^q8(4.29)  and  (4  30)  are  of  the  -  , 

Lj 


same  form  as  in  a  homogeneous  system,  except  for  the  overage  of  recover” 
the  unit  cell. 

terms  of  the  dimensionless  oblique  coordinates  X,  Y, 


(4  30)  becomes  f?>  -  2.\ 


lx  -  (  dX,  j’dY,  r dXj  r  dY^Otj-X,)^  Tj)  -  -  -  (4  3 1  a) 

0  0  -00  -00  T  2itosin^ 


lY-l!«il'dY,  rdX2pdY2(Y2~Y,)^9‘2V,,r2)-  - -  (43^b) 

0  0  -00  -00  T  2nbsin4)  ~ 


fl  ,'l  fOO  fOO  ,  ,  QnCOS^)  . 

IXY  -  j  dX,  j  dY,  J  dXj  j  dYjCXj-XjKY  -Y,)  rj  =  - - (4  3^c) 

0  0  -oo  -oo  T  2Jtstn4i 

In  Appendix  A,  we  show  that  these  sum  rules  (4  31)  ere  Indeed  obeyed  For 

121 

this  purpose,  it  is  convenient  to  express  gy  in  terms  of  the  projector  P  in 
its  form  (49). 

The  above  mentioned  proof  applies  to  any  pehodii^  background, 
including  the  special  case  of  a  lattice  of  fixed  charges  of  negligible  size. 
For  this  limiting  cose,  however,  we  con  also  give  on  alternative  proof, 
which  is  described  in  Appendix  B. 


UoTo 

r  iH  tfWulwie.  we  have  assumed  the  value  1  for  g^A,  the  average  number 
of  mobile  particles  per  wit  cell.  If  g^A  •  p/q  (p,  q  integers),  choosing  a  unit 
cell  q  times  larger  reduces  the  problem  to  g^A  -  p,  end  it  is  easy  to  see  that 
tiMfciflttec  is  solved  by  diagonalizing  a  q  x  q  matrix. 

Whet  about  irrational  values  of  ^A  7  In  terms  of  a  mode)  of  a  crystal, 
it  would  be  a  rather  academic  situation,  since  -^A  is  the  ratio  between  the 
charge  of  an  ion  and  ttw  charge  of  an  electron.  Nevertheless,  this  is  a 
mathematically  interesting  situation.  Furthermore,  our  problem  is  closely 
related  to  a  magnetic  analog^’ of  importance  for  the  theory  of  the 
quantum  Hall  effecr  '  and  also  to  problems  which  arise  in  the  theories  of 
incommensurate  structures^ ' 

We  have  not  been  able  to  compute  the  densities  for  irrational  values  of  I 
g^A.  We  only  went  to  point  out  that,  in  the  simplest  case,  the  problem  I 
reduces  to  studying  the  solutions  of  an  almost-  Mathieu  equation. 

We  start  again  with  the  functions  (4.2),  for  a  square  unit  cell  of  side  a,  | 
and  a  potential  modulation  of  the  form  (4. 13)  . 
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=  [1  +  ‘(cos  2n)(  ♦  cos2n  exp(-ny(X  -  ♦ 


2ni(C*n)Y] 


(4  32) 


where  y  =  Computing  the  projector  P  on  the  space  of  these 

functions  amounts  to  diagonalize  the  matrix  formed  by  the  scalar  products 


Ur  vj  „(r)  -  o(C  -  C  )  A, 


(4.33) 


A  simple  calculation  gives 


Aran  -  (St))"  '"J '  •Ae-x/St  cos  )Sm„ .  ^  Xe-"/2W6„  n.- 1  ‘  6„  , )( 

(434) 


The  problem  of  diagonalizing  A^^^  leads  to 


^Un^  1  ♦  Up. , )  ♦  [cos 


(4.35) 


If  y-1,  the  solution  is  exp(-2niTin),  and  we  retrieve  (4  4)  If  p  is  irrational, 
(4.35)  is  the  almost-Mathieu  equation  in  its  full  glory,  and  we  leave  the 
computation  of  the  projector  P  as  an  open  problem. 


5.  CONCLUSION. 


At  r  -  2,  we  hove  obtained  on  exact  solution  for  the  equilibrium 


statistical  mechanics  of  the  model  of  fixed  ions  and  mobile  electrons 


introduced  by  Hansen  et  or  '  ^  we  have  shown  that,  at  r  -  2.  the  model 


exhibits  the  features  of  a  conducting  phase  ;  the  correlations  at  large 
separations  decay  faster  than  any  Inverse  power  law,  and  the  system  has 
good  screening  properties  (the  Stillinger-Lovett  rule  is  obeyed) 

On  the  basis  of  computer  simulation  results,  it  has  been  claimed  by 


’.-VI  •v\^  ’VI  ’ST. 
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Hansen  and  his  collaborators  that  the  coupling  r^,  at  which  the 
conductor-dielectric  phase  transition  occurs,  go^t6'’^y”u{ji{{eiP vatueTw 
the  radius  of  the  ions  goes  to  zero  Our  exact  results  at  r  =  2  are  not  in 
contradiction  with  this  claim. 


APPENDIX  A  ;  STILMNGER-LOVETT  SUM  RULE 


We  check  the  sum  rules  (4  3 1 ),  where 
r^)  =  -<  r,  IPI  IP|  r,  > 

using  for  <  r,  IPI  >  the  representation  (4  9). 
We  first  consider 


(A  1) 


n' m' 


X  exp{-2ni[Ti(n-m)-T|  (n‘-m  ))-  (X,-C-n)^-  -S-  (x^-^  -n  )^  ♦  2ni(n-n  )Y,} 

X  r* 


oo 

X  I  dX,  exp(-  -IL  (X--C-m)^  -  fl  ()L-C  -m  )^l 

-oo  X*  ^ 

X  1  ”  dY^CYj-Y,)^  exp(-2^(r2)-2nl(C-C  )(Y2-Y,)-2rti(m-m  )Y2]  (A2) 

-oo 


Using  the  periodicity  of  sf,  we  find  for  the  integral  over  Y^  in  (A2) 

■  Tib  ^  ^  'm  .m.M  I'  "''2  «P(-2«(rj>  • 

(2rtr  oc,  0 


In  (A2),  we  can  replace  m'  by  m  ♦  II  end  n'  by  n  ♦  N.  Using  the  periodicity  of 
0,  we  can  replace  X^  -  m  by  X^  and  perform  the  sum  over  m  which  gives  a 

I  oo 

6('n  -  Tl');  we  con  also  replace  X,  -  n  by  X,  and  ^  j  dX,...  by  (  dX,...  We  obtain 

no  -oo 


Iy  - 


eS 

(2n7 


1 

11) 


X  ^  £  exp(2niTiN)  \  dx  expi  -  -5-  (X  -t,-N)^]  |  dY,  exp(-20(r,)-  2niNY,] 

Po  N  Q 


X  -!— £  exp(-2niTiM)  f* dXj  exp(  -  JL  (X^-C-M)^]  /  dY^  ex[i-20irj*  2niMY2] 

Po  f-|  -oo  ^  0 

X  — L_  exp(- J  (X,  -  -  JL  (X- -  0^])  (A3) 


In  (A  3),  as  a  function  of  X,  the  derivative  (3^/3C^)...  is  a  combination  of 
terms  of  the  form  exrt-(n/t)(X,-C)^]  or  (X,-0exp(-(n/x)(X,-0^]  or  (X,-C)‘ 
exp(-(n/x)(X,-O^J  These  terms  enter  the  integral  over  X,  and  generate 


^  I  exp(2niTiN)  Pdx,  exrf-  2  (x,-C)^-  (X,-C-N)^]  (  dY,  exrt-20(r, 

00  N  *°®  ^  t*  0 


-  2«iNY,] 


multiplied  on  the  right  by  1  or  X,  -  C  or  (X,  -  if,  i.e.  either  f(C,  t|)  as  given 
by  (4.0),  or  combinations  of  0f/ac,  af/9Ti,  and  3^f/acaT|.  The  same 
identifications  con  be  mode  with  exp(-(n/x”)(X2-C)^]  and  its  derivatives 
The  result  of  these  identifications  is 


(2nr  0  0  x»x*^ 


I  2L  ♦ 

a-q  (x*xT  f  ail  (x^x*r 


1  21  ]}  (A4) 
f  ac 


Since  f  is  periodic  in  i  and  ii  with  periods  I,  the  contribution  from  the 
derivative  (a/aii)...  vanishes,  and  we  obtain  (4.31b). 

The  computation  of  follows  the  some  lines.  Instead  of  (A3)  we 


obtain 


_V  expCniTiN)  dx,  exp(-  !L(x,-C-N)  ]J  dY,  exp(-20(r,)  -  2ni  NY,] 
^0  -oo  '  r*  0 


X  -1— X  exD(-2itiT]M)  I  dX-  exp(  -  -2-  (X^-C-fl)^]  |  dY,  exp[-20(rj*  2«i 

Po  ri  ^  0 


X^  -  X.  -  ^  {_l_exp(-II(X,  -  0^  --^^(X,  -  0"]}  (A5) 

^  '  2n\  an  ac  fCC.n)  ^  t* 


2n\  ar)  ac,  f(c.n) 


Again,  we  manage  to  recognize  derivatives  of  f(f.,  ^),  with  the  result 


(431c) 


Exchanging  the  X  and  Y  axes  in  (4.3 1  b)  proves  (4.3 1  a) 


APPENDIX  8  ;  LATTICE  OF  POINT  PARTICLES 


We  take  as  the  periodic  background  a  lattice  of  fixed  charged  particles  and 
consider  the  limiting  case  of  point  particles.  We  give  a  direct  proof  that  the 
Stillinger-Lovett  rule  is  obeyed. 

Let  2^  =  Xj  ♦  iy^  be  the  complex  number  which  defines  the  position  of 

the  Ith  mobile  particle;  similarly  let  Zj  define  the  position  of  the  j  th  fixed 

particle.  At  r  »  2.  the  Boltzmann  factor  of  a  system  of  N  mobile  and  N  fixed 
point  particles  is 

Tt  (2|-  2,()  mZj- 

g  ^  ’  J  (B.i) 

where  the  second  form  is  obtained  by  using  an  algebraic  identity  of  Cauchy. 
In  order  to  ovoid  short  distance  divergences,  we  introduce  a  cutoff  at  some 
small  distanced,  and  replace  (B.I)  by 

exp(-pH)  -  l2N|  det  { - 1 - j. ,  nP  (b.2) 

Zj-Zj 

The  limit  a  0  will  be  taken  at  the  end  of  the  calculation.  (B.2)  is  again  a 
squared  determinant,  and  therefore  the  n-body  truncated  densities  are  again 
of  the  form  (2.5),  where  now  P  is  the  projector  on  the  space  spanned  by  the 
functions  (z-Zj)"’{l-exp(-lz-z//2cr^]}.  An  orthogonal  basis  for  this  space 

is  obtained  by  butMtng  Bloch  functions.  Assuming  for  simplicity  that  the 
fixed  particles  are  on  the  square  lattice  Z  -  m  ♦  in,  m,  n  €  Z,  we  define  the 


=  Z  exp(2ni(J,m+'nn)] 
'  mn 


Using  the  Founer  transform 


l-expMCx-mr  ♦  (y-nrl/2o^} 
X  ♦  ly  -(m  ♦ in) 


foo  .00  1 -exp(-(x^*y^)/2cr^] 

g(C,-n)  -  J  dx  J  dy  exp(2ni(C.x>Tiy)]  - 

-oo  -00  x+jy 

exp(-2n^o^(C^Ti^)] 


and  the  Poisson  summation  formula,  we  can  rewrite  ¥  as 


'♦'r  Ti(x,  y)  -  -X  exp{2ni{(C  ♦  m)x  ♦(Ti^n)y]}  g(C^m.  ii^n)  (^5) 
'  mn 


and  we  find  for  the  projector 


<r,  IPI  r,>  -  1  _dU  dt,  -  1  dii 


-1/2  -1/2  I  dxl  dy|¥r  f,(x.  y)r  -i/2  -i/2 

-1/2  -1/2  '  ' 


2^  exp{2ni((C*fn,)x,*(ii*n,)y,-(C*m2)x2-(Ti+n2)y2]}g(C*m,,‘n4n,)g(C*m2' 

,  n,  m^n^ 


Z  ‘n*n)g(C*m.  Ti+n) 

mn  ( 


The  structure  of  (B.6)  allows  us  at  once  to  check  the  overage  density  sum  rule 


f1/2  f1/2  fl/x  rf/iC 

J  dx  j  dy  g(r)  ■  J  dx  J  dy  <  r  IPI  r  >  - 
-1/2  -1/2  -1/2  -1/2 


f1/2  r1/2 


and  the  neutrality  sum  rule 


f 
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1  dX2  J  dyj  fj)  = -1  dXj  f  dyj  <  r,  |P|  Tj  ><  Tj  |P|  r,  > 

-oo  -oo  T  -oo  -oo 


=  -  <r,  |P|  r,  >  = -p(r,)  (B.8) 

We  turn  to  the  Stillinger-Lovett  second  moment 

f  oo  foo  f1/2  fl/2  (2*  _ 

1-2)  dxj)  dy^l  dxj  du,  p^  H  r,,  r2)(x,-X2) 

-OO  -OO  -1/2  -1/2  T 

rOO  fOO  r1/2  fl/2  * 

=-2  J  dx,  J  dyjJ  dx,  J  dy,  <r,  1P|  r^xr^  IPI  r,  >  (Xj-x^) 

-oo  -oo  -1/2  -1/2 

(B.9) 

Using  the  representation  (B.6)  of  the  projector,  we  perform  the  space  integrals  with 
variables  r,  and  Tj  -  r,.  we  find 


\  A/2  rl/2  A/2  fl/2  '  fM;C,tj)f(C.‘n;C.'n) 

I  -  — ^  J  dC  1  dtj  J  dC*  J  dij‘  6*(C-C’)  6(T1  -  Tj  ) 


—2  J  1  . 

2jr  -1/2  -1/2  -1/2  -1/2 


(B.IO) 


when 


(BID 


n.  C.  t|')  -  Z  g(C  ♦  m.  n  ♦  n)  ♦  m,  -n’  ♦  n) 


Replacing  6*(C  -  C  )  by  8(C  -  CK8*/9C‘*),  we  obtain 


1  fi/2  ri/2  I  8*f(C,‘n;C.‘n)  1  3t(M;C'.il) 

I  •  -  —y  J  d(  1  dn ( -  -  ■  - -  ■■■ 

^  -1/2  -1/2  f(w,Ti)  ac8C' 

aftc^Tijc-n) 


h< 

(B  12) 


The  function  f(C,  t);  C,  t|)  defined  by  (B.4)  and  (B.l  I)  is 

expj-4j(*(r^C  ♦  m/  ♦  (i|  ♦  nfn 


t(C,  11,  C.  1})  -  Z - - - 3 — - T 

mn  (C  ♦  mr  ♦  (11  ♦  nr 


(B.l  3) 
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It  IS  convenient  to  display  the  singularity  at  C  =  t)  =  0,  conning  from  the  term 
m  =  n  =  0,  and  to  represent  f,  in  the  integration  domain  lf.|,  |tiI  <  1/2,  by  its 
Laurent  expansion 


T|,  L  T))  =  -2-1-2-  ♦  A  ♦  ... 

c.  ♦Tl 


(B.14) 


In  the  point  particle  limit,  d  -*•  0,  the  sum  in  (B.  13)  diverges  for  large  (m,  n) 
and  A  becomes  infinite,  however  the  terms  of  higher  order  in  (^,  t)) 
represented  by  dots  in  (6.14)  remain  finite  and  can  be  neglected.  If  a  similar 
analysis  is  performed  for  the  derivatives  of  f  appearing  in  (B.12),  it  is  easily 
seen  that  only  the  term  m  «  n  «  0  plays  a  role  as  a  0,  because  the  sums  on 
(m,  n)  which  define  these  derivatives  do  njil  diverge  for  large  (m,  n); 
therefore 


3^U(;,t);Cti)  I  I  3g(^,‘n)  |2 


(B.15) 


af(C,Ti;;’,^)  af(c,'n;C-, 


^'1  .  ~  m..6) 

ac  (,-.1  ^ 


Using  (B.14),  (B.15),  and  (B.16)  in  (B.12),  we  find 


1  f_l/2  f-l/2  A 

I  ”  m  I  I  dtl  '  m  m  'I 
rr  -\/2  -1/2  [WA(C  >1)  )] 


(B.17) 


As  A  goes  to  infirnty,  the  integration  domain  in  (B.17)  can  be  extended  to 
infinity  and  the  integral  evaluated  in  polar  coordinates  Thus 


lim  I--1  (0.10) 

^-*■0  ^ 

This  proves  that  the  Stillinger-Lovett  rule  holds  in  the  point-particle  limit  a 
-*  0 

We  have  not  been  able  to  compute!  for  a  finite  value  of  cr.  It  should  be 
remarked  that  the  Boltzmann  factor  (6.2)  corresponds,  when  d  is  non  zero,  to 
a  complicated  many-body  interaction  which  would  become  the  Coulomb  law 
only  if  all  particles  were  far  apart  from  one  another.  There  is  no  obvious 
reason  for  believing  that  such  a  system  obeys  or  does  not  obey  the 
Stillinger-Lovett  rule.  This  is  in  contrast  with  the  case  of  a  bona  fide 
system  with  extended  fixed  particles,  which  does  obey  the  rule,  os  shown  in 
Appendix  A. 
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Fig.1.  Potenti^modulation  (a),  background  density  (b),  and  particle  density 


(c). 


